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ABSTRACT 

Fluctuations in the 21cm brightness from cosmic hydrogen at redshifts z > 6 were sourced by the primordial 
density perturbations from inflation as well as by the radiation from galaxies. We propose a method to separate 
these components based on the angular dependence of the 21cm fluctuation power spectrum. Peculiar velocities 
increase the power spectrum by a factor of ~ 2 compared to density fluctuations alone, and introduce an angular 
dependence in Fourier space. The resulting angular structure relative to the line of sight facilitates a simple 
separation of the power spectrum into several components, permitting an unambiguous determination of the 
primordial power spectrum of density fluctuations, and of the detailed properties of all astrophysical sources 
of 21cm fluctuations. We also demonstrate that there is no significant information to be gained from 21cm 
measurements on large angular scales. Thus, observational analyses can be confined to small angular areas 
where projection effects are negligible, and theoretical predictions should focus on the three-dimensional power 
spectrum of 21cm fluctuations. 

Subject headings: galaxies: high-redshift, cosmology: theory, galaxies: formation 



1. INTRODUCTION 

Following the recombination of protons and electrons less 
than a million years after the big bang, the universe was 
filled with neutral hydrogen (H I). Hundreds of millions of 
years later, the first galaxies began to reionize the cosmic gas 
(jBarkana & Loeb 2001). The spectra of the farthest quasars 
\ indicate that reionization completed at a redshift z ~ 6, a 
, billion years after t he big bang (Fan et al. 2002; Whit e et all 
■ l2003HWvith e et al.'^2004). The hvperfine spin-flip transition 
' of H I at a wavelength of 21 cm provides the most promis- 
ing tracer of the cosmic gas before the end of reionization. 
Several groups are currently constructing low-frequency ra- 
dio arrays capable of detecting the diffuse 21cm radiation 
(http://space.mit.edu/eor-workshop/). Since this radiation in- 
teracts with H I through a resonant transition, any observed 
wavelength selects a redshift slice of the universe. Future 
measurements of the redshifted 21cm brightness as a func- 
tion of wavelength and direction sho uld provide a three- 
dimensional map of the cosmic H I (Hogan&Reesl lT979t 
iMadau et al. 1997). Fluctuations in the 21cm brightness are 
sourced by primordial density inhomogeneities on all scales 
down to the cosmological Jeans mass, making 21cm the rich- 
est co smological data set on the sky (Loeb & Zaldarriaga 
l2?)0l . 

The 21cm signal can be seen from epochs during which 
the gas was largely neutral and deviated from thermal equi- 
librium with the cosmic microwave background (CMB). The 
signal vanished at redshifts z > 200, when the residual frac- 
tion of free electrons after cosmological recombination kept 
the gas kinetic temperature, T^, close to the CMB tempera- 
ture, r^. But during 200 z -c^ 30 the gas cooled adiabati- 
cally and atomic collisions kept the spin temperature of the 
hyperfine lev el population below r-y, so that the gas appeare d 
in absorption dScott & Reesll990HLoeb & Zaldarriaffal2004 . 
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As the Hubble expansion continued to rarefy the gas, radia- 
tive coupling of Ts to began to dominate and the 21cm 
signal faded. When the first galaxies formed, the UV pho- 
tons they produced between the Lya and Lyman limit wave- 
lengths propagated freely through the universe, redshifted into 
the Lya resonance, and coupled and Tk once again through 
the Wouthuysen-Field (Wouthuysen 1952; Field 1958) effect 
by which the two hyperfine states are mi xed through the ab- 
sorption and re-emission of a Lya photon ( iMadau et al J 19971 
Ciardi & Madau 2003). Emission of UV photons above the 
Lyman limit by the same galaxies initiated the process of 
reionization by creating ionized bubbles in the neutral cos- 
mic gas, while X-ray photons propagated farther and heated 
Tk above throughout the universe. Once grew larger than 
T-y, the gas appeared in 21cm emission. The ionized bubbles 
imprinted a knee in the power spectrum of 21cm fluctuations 
(jZaldarriaga et al. 2004), which traced the H I topology until 
the pro cess of reionization was completed jFurianetto et alJ 
l2004i . 

The various effects that determine the 21cm fluctuations 
can be separated into two classes. On the one hand, the 
density power spectrum probes basic cosmological parame- 
ters and inflationary initial conditions, and can be calculated 
exactly in linear theory. On the other hand, the radiation 
from galaxies, both Lya radiation and ionizing photons, in- 
volves the complex, non-linear physics of galaxy formation 
and star formation. If only the sum of all fluctuations could 
be measured, then it would be difficult to extract the separate 
sources, and in particular, the extraction of the power spec- 
trum would be subject to systematic errors involving the prop- 
erties of galaxies. We show in this Letter that the unique three- 
dimensional properties of 21cm measurements permit a sep- 
aration of these distinct effects. Thus, 21cm fluctuations can 
probe astrophysical (radiative) sources during the epoch of 
first galaxies, while at the same time separately probing phys- 
ical (inflationary) sources. In order to affect this separation 
most easily, it is necessary to measure the three-dimensional 
power spectrum of 21cm fluctuations, while full-sky observa- 
tions are normally restricted to angular projections. However, 
we show that unlike the analogous case of CMB fluctuations, 



2 



21cm measurements on large angular scales are not useful. 
Thus, observational analyses can be confined to small angular 
areas where projection effects are negligible. 

2. SPIN TEMPERATURE HISTORY 

The spin temperature is defined through the ratio be- 
tween the number densities of hydrogen atoms in the excited 
and ground state levels, ni/«o = (^i /^o)exp {-r^/r,} , where 
subscripts 1 and correspond to the excited and ground state 
levels of the 21cm transition, igi/go) = 3 is the ratio of the 
spin degeneracy factors of the levels, and = 0.0682K corre- 
sponds to the energy difference between the levels. As long as 
Ts is smaller than the CMB temperature = 2.725(1 +z) K, 
hydrogen atoms absorb the CMB, whereas if T, > they emit 
excess flux. In general, the resonant 21cm interaction changes 
the brightness temperature of the CMB by ( Scott & Rees 
[1^90; Madau et al. 1997) 7i = r (Z-T-y) /(l+z), where the 
optical depth at a wavelength A = 21cm is 

3cA^Mio«H 

T = ; Xhi , ( 1 ) 

32TTkBTAl+z)(dvr/dr) ' 

where hh is the number density of hydrogen, Aio = 2.85 x 
10"'^ s"' is the spontaneous emission coefficient, xhi is the 
neutral hydrogen fraction, and dvr/dr is the gradient of the 
radial velocity along the line of sight with Vr being the phys- 
ical radial velocity and r the comoving distance; on average 
dvr/dr = H(z)/{1 +z) where H is the Hubble parameter. The 
velocity gradient term arises since the 21cm scattering cross- 
section has a fixed thermal width which translates thro ugh the 
redshift factor ( 1 + v^/c) to a fixed interval in velocity (iSobolevI 
[i960). 

For the concordance set of cosmological parameters 
iSpergel pt al. 200 3). the mean brightness temperature on the 
sky at redshift z is 7), = 28 mK [( 1 + z)/ 1 0] ' [(r^ - T^) /Z] f hi, 
where i'Hi is the mean neutral fraction of hydrogen. The 
spin temperature itself is coupled to T^. through the spin-flip 
transition, which can be excited by collisions or by the ab- 
sorption of Lya photon s. As a resu lt, the combination that 
appears in Th becomes (lFieldlll958l) (Ts-Ty)/Ts = [ji:tot/(l + 
-'^tot)] {^~T^/Tk), where Xtot = Xa+x^ is the sum of the ra- 
diative and collisional threshold parameters. These parame- 
ters aiexa =4PaTi,/27AioT-f and x^ = 4Ki-^)(Tk)nHTj,/3A[()T~f, 
where Pa is the Lya scattering rate which is proportional to 
the Lyg intensity, and K| p is tabulated as a function of 
iAllison & Dalgarno'1969VZvgelman^2Q04V The coupling of 
the spin temperature to the gas temperature becomes substan- 
tial when Xtot 

3. BRIGHTNESS TEMPERATURE FLUCTUATIONS 

Although the mean 21cm emission or absorption is difficult 
to measure due to bright foregrounds, the unique character 
of the fluctuations in Ti, allows for a much easier extraction 
of the signal (Gnedin & Shaver 2003; Zal darriaga et al., 2004; 
[Morales & Hewitt 2004; Morales 2004; ISantos et alJl2004i) . 
In general, the fluctuations in Ti, can be sourced by fluctua- 
tions in gas density, temperature, neutral fraction, radial ve- 
locity gradient, and Lya flux (through Xq). Adopting the no- 
tation 6a for the fractional fluctuation in quantity A (with a 
lone S denoting density perturbations), we find 



Xc dlog{Ki-o) 



where the adiabatic index is 7^ = \ + (8j^lS), and we define 
itot = (1 +^tot)-'i:tot- Taking the Fourier transform, we obtain the 
power spectrum of each quantity; e.g., the total power spec- 
trum Pj,^ is defined by 

(57i(ki)5r.(k2)) = (27r)3<5"(ki +k2)Pr„(ki) , (3) 

where 6j,()&) is the Fourier transform of k is the comoving 
wavevector, and (• • •) denotes an ensemble average. 

4. THE SEPARATION OF POWERS 

The fluctuation Jj,, consists of a number of isotropic sources 
of fluctuations plus the peculiar velocity term -5d,v,- Its 
Fouri er transform i s simply pro portional to that of the density 
field JKaiseil987tlBharadwai & Ali 2004a), 
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where /i = cos 9k in terms of the angle 9^ of k with respect 
to the line of sight. The /i^ dependence in this equation re- 
sults from taking the radial component (oc /i) of the peculiar 
velocity, and then the radial component (oc /x) of its gradient. 
Intuitively, a high-density region possesses a velocity infall 
towards the density peak, implying that a photon must travel 
further from the peak in order to reach a fixed relative redshift, 
compared with the case of pure Hubble expansion. Thus the 
optical depth is always increased by this effect in regions with 
5 > Q. This phenomenon is most properly termed velocity 
compression. 

We therefore write the fluctuation in Fourier space as 

J7;,(k) = /^^^(k) + /3J(k) + 4d(k) , (5) 

where we have defined a coefficient (3 by collecting all terms 
oc (5 in equation (|2j, and have also combined the terms that 
depend on the radiation fields of Lya photons and ionizing 
photons, respectively. We assume that these radiation fields 
produce isotropic power spectra, since the physical processes 
that determine them have no preferred direction in space. The 
total power spectrum is 

PtS^) = iM^Psik) + lAPPsik) + Ps-,M + 

[P'Psik) + P,,a{k) + 2pPs-M] , (6) 
where we have defined the power spectrum P^-rad as the 
Fourier transform of the cross-correlation function, 

6-rad(r)=((5(ri)4ad(ri+r)). (7) 

We note that a similar anisotropy in the power spec- 
trum has been previously derived in a different context, that 
of measuring density fluctuations from galaxy redshift sur- 
veys. In that case, the use of redshifts to estimate distances 
changes the appare nt density of galaxies along the line of 
sight (Kais e3ll987tlL ilie & Efstathiou 1989; Hamilton 199^ 
iFisher. Scharf. & Lahav.19 94). However, in the case of galax- 
ies there is no analog to the method that we demonstrate below 
for separating in 21cm fluctuations the effect of initial condi- 
tions from that of later astrophysical processes. 

We may now calculate the angular power spectrum of the 
brightness temperature on the sky at a given redshift - cor- 
responding to a comoving distance ro along the line of sight. 
The brightness fluctuations can be expanded in spherical har- 
monics with expansion coefficients aim{v), where the angular 
power spectrum 

d^k 



Ci(ro)={\a,Jv)Y)='X^ 



Gf{kro)Psik) + 



Tk-Tj itot dlogiTk) 



(2) 



2Ps-,,d(k)G, (kro)], (kro) + PMjtikn,) 



(8) 
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Ji(x) + (J3 - l)ji(x) and Ji( x) being a linear 



jharadwai 



with Gi(x) 

combi nation of spherical Bessel functions 
l2()04al) . 

The velocity gradient term has been previously neglected 
in 21cm calculations, except for its e ffect on the sky- 
averaged power and on radio visibilities dToz zi et al]l200Ct 
iBharadwai & AM 2004a b). In the simple case where /3 = 1 
and only the density and velocity terms contribute, the veloc- 
ity term increases the total power by a factor of ((1 + /i^)^) = 
1.87 in the spherical average. However, instead of averaging 
the signal, we can use the angular structure of the power spec- 
trum to greatly increase the discriminatory power of 21cm 
observations. We may break up each spherical shell in k 
space into rings of constant fi and construct the observed 
PTi^kjfi). Considering equation (|6} as a polynomial in fi, i.e., 
/x^f^4 + /i^P^2 +P^o, we see that the power at just three values 
of /Lt is required in order to separate out the coefficients of 1 , 
fi^, and fi'^ for each k. 

If the velocity compression were not present, then only the 
/i-independent term (times T^) would be observed, and its 
separation into the five components (r^, and three power 
spectra) would be difficult and subject to degeneracies. Once 
the power has been separated into three parts, however, the /i'' 
coefficient can be used to measure the density power spectrum 
directly, with no interference from any other source of fluctu- 
ations. Since the overall amplitude of the power spectrum, 
and its scaling with redshift, are well determined from the 
combination of the CMB temperature fluctuations and galaxy 
surveys, the amplitude of directly determines the mean 
brightness temperature Th on the sky, which measures a com- 
bination of Ts and xui at the observed redshift. Once Ps{k) 
has been determined, the coefficients of the fi^ term and the 
/^-independent term must be used to determine the remain- 
ing unknowns, f3, Ps-ndik), and Prad(^)- Since the coeffi- 
cient f3 is independent of k, determining it and thus break- 
ing the last remaining degeneracy requires only a weak addi- 
tional assumption on the behavior of the power spectra, such 
as their asymptotic behavior at large or small scales. If the 
measurements cover Nk values of wavenumber k, then one 
wishes to determine 2Nk + 1 quantities based on IN^ measure- 
ments, which should not cause significant degeneracies when 
Nk ^ I- Even without knowing f3, one can probe whether 
some sources of Prad(^) are uncorrected with 6; the quantity 
P,,.s(k) = P,,o-p2^/(4P^4) equals P,,d-Pl,JPs, which re- 
ceives no contribution from any source that is a linear func- 
tional of the density distribution. 

5. SPECIFIC EPOCHS 

At z ^ 35, galaxy formation is in its infancy, while 
collisions are effective due to the high gas density 
(p..oeb & Zaldarriaga 2004). Thus, Xtot = Xc ^ Xa, xhi = 1, 
and the gas cools adiabatically (with ja = 5/3). In addition 
to me asuring the density power spectrum dLoeb & Zaldarriaeal 
120041) . one can also trace the redshift evolution of «hi, T-^, and 
Tk, and verify that collisions are causing the fluctuations. 

At z < 35, collisions become ineffective but the first stars 
produce a cosmic background of Lya photons that couples T, 
to Tk- During the period of initial Lya coupling, fluctuations 
in the Lya flux translate into fluctuations in the 21cm bright- 
ness jBarkana & Loebl l2004'). This signal can be observed 
fromz ~ 25 until the Lya coupling is completed (i.e., jctot ^ 1) 
at z ~ 15. At a given redshift, each atom sees Lya photons that 
were originally emitted at earlier times at rest-frame wave- 
lengths between Lya and the Lyman limit. Distant sources 
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Fig. 1 . — Observable power spectra during the period of initial Lya cou- 
pling. Upper panel: Assumes adiabatic cooling. Lower panel: Assumes 
pre-heating to 500 K by X-ray sources. We show P^4 = Pg (solid curves), 
P^2 (short-dashed curves), and Pun-^ (long-dashed curves). We also show for 
comparison 2/3Pg (dotted curves). 



are time retarded, and since there are fewer galaxies in the 
distant, earlier universe, each atom sees sources only out to an 
apparent source horizon of ~ 100 comoving Mpc at z 20. 
A significant portion of the flux comes from nearby sources, 
because of the l/r^ decline of flux with distance, and since 
higher Lyman series photons, which are degraded to Lya pho- 
tons through scattering, can only be seen from a small redshift 
interval that corresponds to the wavelength interval between 
two consecutive atomic levels. 

There are two separate sources of fluctuations in the 
Lya flux ( B arkana & Loeb 2004) . The first is density inho- 
mogeneities. Since gravitational instability proceeds faster in 
overdense regions, the biased distribution of rare galactic ha- 
los fluctuates much more than the global dark matter density. 
When the number of sources seen by each atom is relatively 
small, Poisson fluctuations provide a second source of fluc- 
tuations. Unlike typical Poisson noise, these fluctuations are 
correlated between gas elements at different places, since two 
nearby elements see many of the same sources. Assuming 
a scale-invariant spectrum of primordial density fluctuations, 
and that Xq = 1 is produced at z = 20 by galaxies in dark matter 
halos where the gas cools efficiently via atomic cooling, we 
show in Figure 1 the predicted observable power spectra. The 
figure shows that (3 can be measured from the ratio P^i /P^t 
at A; > 1 Mpc"', allowing the density-induced fluctuations in 
flux to be extracted from P^i, while only the Poisson fluctu- 
ations contribute to Pun-<5- Each of these components probes 
the number density of galaxies through its magnitude, and the 
distribution of source distances through its shape. Measure- 
ments at A: > 100 Mpc"' can independently probe Tk because 
of the smoothing effects of the gas pressure and the thermal 
width of the 21cm line. 

After Lya coupling and X-ray heating are both completed, 
reionization continues. Since [3=1 and Tk ^ T^, the nor- 
malization of Pfji directly measures the mean neutral hydro- 
gen fraction, and one can separately probe the density fluc- 
tuations, the neutral hydrogen fluctuations, and their cross- 
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Fig. 2. — Effect of large-scale power on the angular power spectrum of 
21cm anisotropies on the sky. This example shows the power from density 
fluctuations and velocity compression, assuming a hot IGM at z = 12 with 
Ts = Tii ^ Tj. We show the % change in C; if we were to cut off the power 
spectrum above l/ztof 200, 180, 160, 140, 120, and 100 Mpc (top to bottom). 
Also shown for comparison is the cosmic variance for averaging in bands of 
Al I (dashed lines). 



correlation. 

6. FLUCTUATIONS ON LARGE ANGULAR SCALES 

Full-sky observations must normally be analyzed with an 
angular and radial transform, rather than a Fourier transform 
which is simpler and yields more directly the underlying 3D 
power spectrum. In an angular transform on the sky, an an- 
gle of 6 radians translates to a spherical multipole I 3.5/9. 
For measurements on a screen at a comoving distance ro, 
a multipole / normally measures 3D power on a scale of 
r' ~ Oro - 35 // Gpc for Z > 1, since ro ^ 10 Gpc at z > 10. 
This estimate fails at / < 100, however, when we consider the 
sources of 21cm fluctuations. The angular projection implied 



in Ci involves a weighted average [equation (|8}] which favors 
large scales when / is small, but density fluctuations possess 
little large-scale power, and the C/ are dominated by power 
around the peak of kPs{k), at a few tens of comoving Mpc. 

Figure 2 shows that for density and velocity fluctuations, 
even the / = 1 multipole is affected by power at > 200 Mpc 
only at the 2% level. Due to the small number of large angu- 
lar modes available on the sky, the expectation value of C/ 
cannot be measured precisely at small I. Figure 2 shows that 
this precludes new information from being obtained on scales 
A:"' > 130 Mpc using angular structure at any given redshift. 
An optimal observing strategy would be to divide the sky into 
separate angular pixels of size 9 < 2°, in which the curvature 
changes distances by only 6*^/24 ~ 7 x 10"^. This would al- 
low a direct 3D Fourier transform, which is also the most nat- 
ural way of analyzin g radio observations (.Morales & Hewitt 
l2()04llMoralesl2004l) . 

7. SUMMARY 

We have argued that observers should try to measure the 
anisotropic form of the power spectrum of 21cm fluctuations. 
Since measurements on large angular scales are not worth- 
while, the three-dimensional power spectrum may indeed be 
measured by dividing the sky into small angular areas. The 
anisotropy itself, which is due to the effect of peculiar veloci- 
ties, facilitates a simple separation of the power spectrum into 
three separate components. One component allows an unam- 
biguous determination of the primordial power spectrum of 
gas density fluctuations, which if measured will probe basic 
cosmological parameters as well as inflation. The other two 
components probe the properties of astrophysical (radiative) 
sources of 21cm fluctuations. 
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